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TRANSMISSIpN SYSTEM 

(57) Abstract i 



A communications system employing the Discrete Multitone (DMT) line 
code divides the used frequency iband XV into a large number of subchannels, 
each with a bandwidth Af=W/N. ^he signal ih each subchannel is independently 
modulated at the symbol rate 1/A/ Because of noise and the non-ideal 
characteristics of the transmission channel, ; the signal-to-noise ratio (SNR) 
varies for different subchannels across the frequency band. Each subchannel 
is used to transmit an integer number of Sits in a two-dimensional signal 
constellation, the susceptibility : to noise increasing with the size of the 
constellation assuming that the same signal'energy is used. In order to use 
the transmission channel optimally, subchannels are assigned different numbers 
of bits depending on their specific SNK values, i.e., a subchannel with a high 
SNR value is assigned a large number of bits and vice versa. An optimal bit 
allocation is determined that ensures a certain bit error rate (BER) is obtained 
for all subchannels. In order to simplify numerical computations, thresholding 
is employed to make initial bit assignments.' Further simplification of the bit 
loading procedures is obtained by ending bit allocation optimization efforts after 
a targeted BER is achieved. .. 
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WO 99/16224 PCT/SE98/01598 



REDUCED COMPLEXITY BIT ALLOCATION TO SUBCHANNELS IN A 
MULTI-CARRIER HIGH SPEED DATA TRANSMISSION SYSTEM 



FIELD OF THE INVENTION 

The present invention relates to multi-carrier, high speed data 
5 transmission and reception, and more particularly, to a method and system for simply 
and cost effectively allocating optimal bandwidth to each of the multi-carriers. 

BACKGROUND OF THE INVENTION 

The general principle of transmitting data by dividing it into several 
interleaved bit streams, and using these bit streams to modulate several carriers is well 
10 known. Modems using this principle employ what is referred to as Discrete Multitone 
(DMT) modulation and demodulation. These types of modems are being used or 
considered for use in such applications as cellular radio and High rate Digital 
Subscriber Lines (HDSL). 

In a L Discrete Multitorie system, the! input bit stream is first serial-to- 
15 parallel converted. The parallel outpiut is then grouped into N groups of bits 

corresponding to the number of bits^per symbol. Portions of bits are allocated to each 

to the. number of bits transmitted oVjer each DMT subchannel. 

■•J . ..... . ■ 

Like a single carrier system, a multicarrier data transmission system tries 
20 * to maximize a data transmission or bit rate for a particular Bit Error Rate (BER) and a 
particular Signal to Noise Ratio (SNR). In theory, the total Bit rate (transmitted over. all 
of the sub-channels) is maximized if the bit error rates for all the subchannels are equal. 
The concern is how best to allocate bits to each subchannel to attain the same BER over 
all of the subchannels. - ! 
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! Fig. 1 shows an example Discrete Multitone (DMT) communication 
system^ in- which the present invention may be advantageously employed. Transmitter 
10 includes a seriaMo-pfarallel convert6r44 5 a-multicarrier modulator: 16, and a 
pretransmit processor 18. Receiver 12 includes apost channel^processdr 20, a 
5 multicarrier demodulator 22, and a parallel-to-serial converter 24. The transmitter and 
receiver are linked in this example by a digital subscriber -line (DSL) or other form of 
communication channel 26. Serial input data at a rate of btotal^T bits per second are 
grouped by converter 14 into blocks of btotil'bits for each multicarrier symbol, with a 
symbol period of T. The b tota i bits in each multicarrier symbol are used to modulate N 
.10 . , separate carriers in modulator 16 s wtfrbj, bits ? mpdulating ther* carrier. 

A preferred example embodiment uses an Inverse Discrete Fourier 
Transform (IDFT) during modulation to generate N s time-domain samples of a transmit 
s ignal for each block of btotai 'Bits, where ^ s is preferably equal to 2N. The 
corresponding multicarrier demodulator performs a Discrete Fourier Transform (DFT), 
15 where b; bits are recovered from the i"" 1 carrier. As depicted in Fig. 2/the carriers or 
subchannels in a DMt system &re spaced 1/t Hz apart across N/T Hz of the "frequency 
band. More detailed discussion of the principles of multicarrier transmission and 
f , reception in general is. given by J.. A. M C.MBmgham,in "Multicagier Modulation for Data 
v .. Transmission; An Idea Whose Time, Has Come", IEEE Communications. Magazine, 
2p .. Volume.28, Number 5, pp. 5-14, May ,1990. . - r ; . 

The general structure of a Discrete Multitone Modulation DMT system is 
Illustrated in Fig. 3, where {^Xj 1 ,...; X^f-l} are the original, complex, input data 
symbols, {x^} is the modulated data sequence (before cyclic prefix)', {h^} is the 
discrete-time channel response, {njj is the additive noise sequence, {y^} is the 
25 received sequence (after the removal of cyclic prefix), and { x o, x \ , x N-l } are the 
decoded, complex data symbols. The pj's and pj* *s in Fig. 3 are known as the 
modulating and the demodulating vectors, and preferably they are chosen to be 
orthonormal. 
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. Adding the cyclic prefix is. a discrete r time technique used to eliminate 
interblock interference (IBI) in the DM! system. - The independent- modulating and 
demodulating vectors in Fig. 3 are the IE>FX»and the DFT vectors, giyen by the 
following pair of relationships: . , r : . ■ 

1 2N— 1 -J2mim 

5 - X ^l2Nnh X " e - • -'V' 0) ' 

■x.—±± l I. X„e - " (2) 

^ : A DMT'System witfr^ in the frequency domain 

requires a DFT size of 2N, and the forced conjugate symmetry in the frequency domain 
results in the desired real-valued time domain samples. In the preferred embodiment, 
10 IDFJ and DFT are implemented with the; well known IFFT and FFT algorithms using 
the. approach disclosed in. commonly-assigned "Method and Apparatus for Efficient 
. Computation of Discrete Fourier Transform and Inverse Discrete Fourier Transform" 
U.S. patent .application Serial No. 0$/887,467 ? (attorney docket 1410-340), filed on July 
2, 1997, t ( ' f ..... ... J.-, , • , t . . , . , , 

15 • 1 " * '•' Fig. 4 is -a more 'detailed 'block diagram showing the principal operative 
- components of a basic DMT transmitter 3 G and a basic DMT receiver 32 connected 
through a channel 34. Serial input datia^re grouped into blocks, converted to a parallel 
form, and appropriately encoded. by an encoder 36. Each one of the N active 
subchannels contains a number of bk bits allocated to that subchannel. The total 

20 number Qf bits in a DMT symbol is therefore: 

. B total = S V . . V (3) 
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,r -> i > : Each of the N bit groups is. mapped by the symbol into a two-dimensional 

-signal constellation. The relation between .the .number, of bits, in each block of bits bfc 

'and the necessary number of signal points*, in- a:- constellation M. is given by: 

: "-< ' c ; - : o j-.^j f-:.-'i.;-j.;;*\* r i . . 

... , , . = 2 6 ' . ' (4) 

The signal points in a ctinstellatioti are typically arranged in a rectangular 
pattern for' even values of bfc and a crbss'fcr odd Allies of % ; Provided that the same 
signal energy is used for all constelfati&isy thedistahdfc betweeri : neighbdring signal 
points is reduced .with increasing consteljation size resulting in increased BER with the 
. .same sjgnal-tornoise ration (SNR). , ; v . t : . v , mU , , 

" ' j The ! butpurfrom the encodfef is complex numbers, one for each group of 
: bits, which are their Ted' to a device that calculaies the Inverse Discrete Fourier 
Transform (IDFT). the output is^teaPs^^^ be considered the '"' 

superposing of ;V modulated <^6^6ii^Gim&cs spaced A/apart/ s 

. . , The parallel IFF J output is c : on^.prted back to a ,seiial data stream by a 

converter 40, The.digital.modulated (jata^tream is cy.clically prefixed, converted to 
analog form by a digital-to-analog converter (DAC), 42, : low-pa$s filtered at 44, and 
passed through a D.C. isolating transformer 46 during pre-transmit processing to 
produce ah* analog signal transmitted over a ti*ahsmikion Channel 34. 

rrr At the receiver end, the received analog signal is passed through a D.C. 
isolating ; transformer; ^nd.low-pass filter 48,,converted to digital form by .an analog-to- 
digital converter (ADC) 50, time domain pre-equalized by a finite impulse response 
(FIR) Jilter 52 to limit the effective memory of the channel, and stripped of the, cyclic 
prefix during post-receive processing in converter 54. The resulting.digital signals are 
demodulated by, 2N FFT operation 56, and .converted to parallel frequency domain 
signals. Since th<? amplitude vs. frequency and the..d?lay vs. frequency responses of the 
channel are not necessarily constant across the entire used frequency band, the received 
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signal will differ from the transmitted signal, and -the parallel inputs to the decoder 58 
will differ from those parallel outputs -from, the encoder 36. One form of equalization 
used to compensate these differencesissa' frequency dofnain equalizer (FEQ) which 
individually adjusts for the attenuation and delay of each of the carriers immediately 
before the parallel frequency domain signals are passed to the decoder. The frequency 
.domain equalized signals a.re appropriately decoded#nd converted back to a serial form 
e by.the decoder 58. Ideally* .the. jdiet^Gted r QUtput serial. data from the decoder 58 is the 
same asthe input serial 4ata to the, encoder r 3 6., . , « . . 

In a conventional' sirigfe-cMw . 
Quadrature Amplitude Modulation (QAM) System implemented with ari equalizer or a 
encoder,, the transn^ssion,bandwidth.is determined by, the symbol rate and the carrier 
frequency of the system^ ynfortunately.,. variable symbol rate, single-carrier systems are 
still, impractical to implement from„a,complexity 5 standpoint with today '$ technology, 
and even if they are implemented, the granularity of possible symbol rates are typically 
very coarse. This is a direct consequence of the fact that for a fixed data rate, theu 
symbol rate can ofily change iri tiisfcrete multip where 
b'symbol is the ritaibero'f bits transmitted by each data symbol, prbvidedlhat only 
signer constellatittns with integer numbers 6f Tbits are used/ 4 

: Multicarrier modulation, however, offers much more flexibility in 

CTrir: - :!! ITl*- * *~ ~ ~ ~ o ~~+~ . "Ul ~ U ~ 1 ~ 4-U ~+ ^ n+ /->-T o mn/'Vi 1 nrrr^y rmmUa* ^-P 

bits jper symbol over a largfe number of cafriers. Since different numbers of bits can be 
' transmitted through the different carriers (subchannels), the multicarrier transceiver can 
contrdl of the transmission bandwidth. As a result, more data bits are transmitted 
through the better carriers/i.e;, having a greater Signal to Noise Ratio (SNR), and less 
data bits are transmitted through the worse carriers, i.e., having a l6\ver SNR as shown 
"in Fig. 5. In this manner, optimal performance can be achieved,- in terms of either 
maximizing total data throughput for a fixed sysfemi peirfoimahce" margin or maximizing 
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' over air system performance margin for a fixed target dat'a rate under a fixed- total input 
• power constraint and. target bit-error-rate,. 

Various bandwidth/bit allocation optimization techniques for multicarrier 
modulation have been proposed for different applications. The theoretical optimal 
energy allocation is the well known water-pouring energy distribution as described for 
example in J,G: Proakis, (1995); Digital Communications 3 rd edition, McGraw-Hill, 
. New York. : Conceptually, the water'pourinf energy distribution model is obtained by 
first-inverting the channel SNR curve 'shown in Fig'. S, in: the frequeneydomain and then 
pouring the available, energy into the- re'sulting- curve [ ; The- final '"water level 5 ' 'at each 
particular frequency/carrier/subchannel;i$' the optimal, amount 6f energy to us£ at that 
particular frequency/carrier/subchannel. The number of bits transmitted by each carrier 
is determined by the maximum number of bits that can be supported by each carrier, 
given the water-pouring energy distribution and the desired bit-error-rate of each 
carrier. While the water-pouring technique yields the optimal bit and energy 
allocations, and therefore the optimal transmission bandwidth, it is typically complex 

and impractical to implement. 

• j! ^ :vV ■,■:/■.;-<■ r OA r v.". 0 ■■ r 

r • o , •'. Another. multicarrier bandwidth optimization technique designed for 
voice-band modem applications, where only constellations 'with an integer- nufnber of 
bits are used, was proposed by D. Hughes-Hartogs and described in "Ensemble Modem 
Structure for Imperfect Transmission Media", European Patent Number 0 224 556 Bl 
(published October 1991). The basic idea behind Hughes-Hartogs' algorithm is to 
incrementally add, one bit at a time, the amount of data to be transmitted in each 
multicarrier symboluntil a desired data fate or power constraint is satisfied. 
Furthermore, this algorithm always chooses the subchannel that requires the least 
amount of incremental ppwer to suppqrt : an additional bit at Jhe given desired bit-error- 
rate. Hughes Hartogs' algorithm, however, js typically too slow in convergence, 
making it especially .unsuitable for applications where there are large numbers of 
carriers atjd bits to be supported by eaph niultiQaxrief symbol, as ppinte^i opt in P.S. 
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Chow, J.M. Cioffi, J.A.C., Bingham, "A Practical! Discrete Multitone Transceiver 
Loading Algorithm for Data Transmission Over Spectrally Shaped; Channels," IEEE 
Transactions on Communications, Vol. 43, No. 2/3/4, pp. 773 775, 1995. Also, the 
channel response is assumed to be fixed, i.e., time-invariant after training, which is not 
a good assumption for many applications, e.g., mobile radio communications. 

- :. In. order to achieve faster convergence, at Jhe expense of a slight • 

performance degradation, Ghowje&ak proposed a three-step bit loading algorithm. 
First, an optimal system ^ performance-margin is found; second; convergence is 
.guaranteed with a suboptimaljoop?: and third, the. energy distribution : is adjusted for 
each subchannel if the transmission .systemr is constrained by the total transmit power. 

However, Chow's bitioading algorithm requires computationally complex 
steps best run on a microprocessor and a number of iterative elements that may be;: ; 
excluded. Both requirements increase the system cost and reduce speed. The bit 
loading procedure in accordance with the present invention minimizes computational 
complexity and permits execution by specific hardware which is faster than a 
programmed microprocessor driver solutions. Also those parts of the algorithm that are 
best run on a microprocessor .are oif. such basic character that only a rudimentary, 
, inexpensive device is necessary; v; - _ ' • ; '. I 

It is an object of this invention therefore to provide optimal bitioading 
using less computational power and at a fast rate. 

SUMMARY OF THE INVENTION 

" * ; Bits are allocated among subchannels in a multi-carrier modulation 

environment, under the constraint of a targeted bit error rate (BER) as well as a targeted 
■-number of bits per Symbol. Thresholding is used to make initial bit assignments to 
subchannels in ordfcr to simplify numerical calculations! If tfre targeted BER cannot be 
achieved with the' targeted tiuxrfber of bits/ the overall 'channel is insufficient, and a 
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higher BER or a lower transmission rate must be considered. To further increase the 
efficiency of the process, a "block-redufction" step may be employed which reduces by 

one the number of bits allocated to each active subchannel if the total number of 

1 - ' : . ■ . • " ' .' i .j t-/ 1 • - -ir, . • fo;; ., -\. - . fc . 

allocated bits exceeds the number of targeted bits by more than the number of active 
subchannels. Further fine tuning may also be performed. 

- : The simplified bit allocation procedures in accordance with the present 

invention" may be performed quite rapidly ^to permit adaptative allocation Of bits in real 
time td accommodate, for 'exampl^j-ehan^irig-cHarinelxohditions. In addition, the 
simplificd'procedures permit irnpleiffentatiOil using dedicated hardware for improved 

'speed':' ■ ; ' ! " -I \Z ':V» ":<..■ ' ■ •■ :.-< j 

:i ; - -Mo're specifically;; in af multi-carrier commuhications system with a 
•transmitter and a receiver linked together by a communications channel affdeted by 
; noise and other channel impairments, "datafcitS iolie transmitted over the channel at a 
■prescribed data transmission rate'and'Mt ferror'ratd are -allocated to plural subchannels. 
A series 6f thresholds is calculated for eacft subchainnel -with each threshold * 
- corresponding to a preset number of bits to be allocated to that subchannel; Each 
subchannel is associated with a signal constellation where bits allocated to 5 the 
subchannel are mapped to signal points in the signal constellation. The subchannel 
thresholds' are calculatedibased oh ; tMprescribed error rate, an average energy and size 
of the subchannel signal constellation. The threshold calculations need only be 
performed once, preferably before the transmitter/receiver is brought on-line. 

Subsequently, after the transceiver is brought into operation, an error 
estimate such as a noise standard deviation is obtained for each subchannel. That error 
estimate is compared to the thresholds calculated" for the subchannel. The number of 
bits is then allocated to each subchannel corresponding^ one of the thresholds 
exceeded by the obtained error measurement. The error estimate- comparing and 
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obtaining procedures as well as tjie bit allocation may be performed in, real time 
whenever the; channel condition changes. , f ,y. ■ . -s • 

. t „ , r. . • .... „ , * 

' A determination's made whether a total number of bits allocated to all the 
subchannels is greater than or equaVto a total target number of bits associated with the 
prescribed data transmission rate. If not, a decision is made whether to reduce the 
prescribed data transmissiqn rate .or to reduce the prescribed bit error rate. Otherwise, 
, the number of bits allocated,to qnje pr^mo reduced until the 

. total .number of allocated bit? is >; n?arly r ;equai to ; the total target number of bits. In 
particular, the number of bits, to hej aUocated to each active subchannel rnay.be reduced 
by one until the excess number of bits is less than the number of active subchannels. 

/},. -, * When ^jQTj^ain^Qr'OfithQ ^.excess nurnber.,ofbit^ is less than the nunjber of 
* active subchannel^, the reniainder may be further reduced by reducing the number of 
bits allocated: to selected ones the, active subchannels. Bits are subtracted based on a 
magnitude of a difference, between the^ error estimate and- the exceeded thresholdyfor 
each active subchannel.. Those actiy^ subchannels having larger difference^magnitudes 
are selected first, far the ; bit reduction followed by subchannels having smaller ; - 
: difference t magn-itudes. v ,. l -,; v ;. :tif - -x ; ■ 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram illustr^tiijg a basic prior art multicanrier 
transmitter, communication channel and receiver; 

. Fig. 2 is a diagram conceptually illustrating frequency domain discrete 
multitone modulation (DMT); K ........ 

'■ Frg. 3 is a diagram schematically depicting the basic structure of a 
1 discrete ^multi tone' modulation system; 
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*r.- Fig. 4*is a.block. .diagraift depicting the' principal components of a discrete 
: niuftitone-.(DMT) modulation communication system; 

Fig. 5 illustrates DMT bit allocation in the frequency domain; 

Fig. 6 illustrates two-dimensional 4 and 16 QAM signal constellations; 

Fig. 7 shows one example configuration implementing the thresholding 
operation and bit allocation procedures in accordance with one example embodiment of 
the present invention; 

• . • Fig. 8 is a .diagram providing £ specific ^example threshold and bit 
allocation operation based on the example configuration shown in Fig. 7; 

' " ; Fig: 9 show's ari alternative example configuration for implementing the 
. thresholding and bit allpcatipn operations in accordance \yith the present invention; and 

Fig. 10 illustrates bit loadihg pr8c£cfures in flowchart form m accordance 
with the present invention. i: '- ,J - 1 - - ' :,Vi| ■-'^ * - ' tl " 

DETAILED DESCRIPTION OF THE INVENTION 

In 'the following description, for purposes of 'explanation and not 
• limitation, specific details are set forth, .such as .particular hardware, applications, 
techniques, etc. in. order to provide a thorough understanding of the present invention. 
However, it will be apparent ; to one.skilled in the art that.the present invention may be 
practiced in other embodiments that depart from these specific details. In other 
instances, detailed descriptions of well known methods, protocols, devices, and circuits 
are omitted so as not to obscure the description of the present invention with 
unnecessary detail. 

The example apparatus that may be used in implementing the present 
invention may include the transmitter and receiver represented by the block diagram of 
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FIG. 4 already described above that in addition includes hardware and/or software for 
implementing the novel bit allocation procedures according to the present invention. 
The theoretical foundation for those bit allocation procedures now follows. 

Theoretical Foundation 

The bitloading procedures in accordance with the invention allocates bits 
to be transmitted in such a way that each subchannel satisfies the demands on 
probability of a symbol error, Pj^, within a system performance margin. This bit 
allocation is based on a noise variance, <r , in each subchannel and prior knowledge of 
the average energy in eidh subchannel Signal constellation, S dv , which together 
effective^ provide the sigrial-to-rioise ratio SNR for eacli subchannel. Any additional 
safeguard, margin of bits is distributed evenly , among the subphannels. 

^ According to j.Gr. PMakis; (1095), Digital Commuriicatibhs 3 r <* edition, 
McGraw-Hill, r New YorJc 5l the symbol ; error probability for a rectangular 
constellation in Additive White Gaussian Noise (AWGN) is tightly upper bounded by: 

where a number ;of constellation-points Mper subchannel is defined as 
M'=2 h ; b = lv2 5 .... 9 6 max r is the number of bits conveyed per each signal point in the 
Constellation; £ av is the average energy of the subchannel constellation anil is a known 
design parameter set for each subchannel constellation; and N 0 is the single-sided 
power spectral density of the nofee for a received signal' constellation point (at the 
receiver) and is defined as A^ 0 = Icr , where cr is the received noise variance per 
real/imaginary dimension. The noise variance can be visualized (in a simplified 
manner) as the variance of the symbol error, i.e., between the signal point received and 
the .closest consteUatipn. signal point in the complex plane. - 
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; ; 1 The Q-fiihction in equation (5) relates the signal-to-noise ratio (roughly 
defined as { e av f-N 0 -) t6 thfesyniboP error Pjvl for a particular constellation size and is 
: ' defined -mathematically as: ' *' '•' - - 1 J- '* :A •" 

Q(x) = ~ r e^du (6) 

where u is the integral variable. In other words, the Q-function of a certain value x is 
calculated^ solving .. .. . 

■ 1 ' Even though Pjvfand 3 BER :i ai'e riot identical, this analysis assumes that the 

symbol error probability arid bit ehxff rdtb : (Bfi£) 'are approximately the skme because 
one symbol usually corresponds to one or more "bits: Certain symbol encoding 
techniques, (i.e.. Gray coding where adjacent symbols, containing the most likely errors 
only differ by one bit), may be employed to ensure that the bit error probability and the 
symbol error probability are very similar. 

■ - ' ' Jn most practical applications;- *he bit error rate (and in this 'application the 
symbol error rate Pm) is a known system design-constraint. To achieve the same BER 
(or.Pjyi) for all subchannel signal constellations and-ther.eby optimally allocate bits to 
each subchannel under a constraint that the average energy of all subchannel 
constellations is the same, the denominator of equation (5); (M-l) No must be 
maintained constant, i.e., ' ' 

. (M-l)^ 0 = 2a 2 (2 A -.l) = C ,(7>, . . .. 

The constant C depends oh the 'specific BER to be achieved and may be calculated 
according to: '■ 
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Solving equation (7) for the number of bits to be allocated, 6, provides the 
number of bits that may be assigned tp a-subchanne] given ascertain QER^and the noise 
variance a] on that particular subchannel, i.e., the bit allocation for each subchannel k 
corresponding to b k : 

where a k is the noise variance perreaF(i) &nd imaginary (Q) dimensioh on 
subchannel k In a preferred implementation,- .the. bit loading values are limited to 
integer constellations. e A g iV by rounding. them Jo integer constellation sizes which affects 
the total number of allocated _l?ijs.. : nf; - ■ 

Iri equation (7). both 6 and & may varyV The subchannel noise variance 
o* A varies because of changing noise on the subchannel. Therefore, the number of bits 
bk allocated to that subchannel must be varied to compensate for the subchannel noise 
phanges.to ensure thatequation (7); is satisfied; as a;constant and equal to C. If the* 
-subchannel- noise- variance :<t;: increases',;^ must 'decrease, and fewer bits are sent'over 
:the subchannel If decreases, cah be increased With more bits being sent. 

As already described, for a particular system design, the parameters 
or BER, Q, and s m , are preset/known. Therefore, only the noise variance a 2 needs to 
be determined for each subchannel to solve for the optimal bit allocation per subchannel 
according to equation (9). To obtain the'subchahnel noise variance a 2 , the noise 
standard deviation a k is determined based on a series of error measurements for 
symbols received on the subchannel. The subchannel noise standard deviation a k is the 
input to the bit loading procedures employed in the present invention and corresponds 
to the mean of the received symbol error magnitude (attributed to noise) for each 
subchannel. 
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- The foll'oWihg described how the noise standard deviation estimate a k 
.. may be determined for each subchannel k.. The received symbol error is the vector 
difference between a received signaj. point and the closest signal constellation point for 
.a particular subchannel. Thus, the received symbol error is a complex number. To 
estimate ^deviation of the error, either the real (I) or the imaginary (Q) part of the 
received symbol error may be estimated for a theoretical ideal channel with truly 
Additive White Gaussian Noise (AWGN) because the I and Q components would have 
_ equal variance, i.e., the actually received signal points would be equally distributed in a 
circular .pattern around the signal constellation point. In. reality, the noise variance may 
be different in the I and ; Q dimensions,, i.e., not a circular.pattern, resulting in poor noise * 
deviation estimates if noise in only one of the two dimensions is measured. 

Thus, it is preferred that the noise magnitudes in both the I and Q 
dimensions be estimated. The received symbol error magnitude can be expressed as: 

. , r . f ■ u v ■-. L>* = 4 X h+A- - v - * > -.(JO) 

■ 4 '•-/. • . ■■•'/.■>; . , ... j ; ih''' 1 ; ) ■ ■ is. r »; - ";v.-'- 
where ^, and X 2 are zero-mean, statistically-independent, gaussian-random variables, 

each having a variance cr 2 corresponding to the noise components in the complex I- and 
Q-dimensions. According to J.G. Proakis, (1995), Digital Communications 3 rd edition, 
McGraw-Hill, New York, the random variable R is Rayleigh distributed with an 
expected (E) mean value: 

<? = £(J?) = ^cr (11) 

From equation ^! 1), the measured noise .standard deviation estimate a k for each 
subchannel is preferably corrected using E(R), (i.e., by dividing j the measured a k by 
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Example Embodiment of Subchannel Bit Allocation Procedures 

The subchannel bit allocation procedures in abcof dance With' th [ e present 
invention take as an input vector "tfe coirecfed noise standard deviation estimate a k for 
each subchannel k. As described above,' the main goal in optimally allocating bits is to 
maintain a constant, predetermined symbol' error probability Pf^/BER across each of the 
subchannels. Towards that goal, tlie invention uses thresholding to make initial bit 
assignments rather than using conventional complicated numerical computations. 
Further simplification is achieved by ending bit aflocation optimization efforts after a 
target P^/BER is obtained. In another example embodiment^ the preferred bit 
allocation procedures includes three b>asic steps that each converge within a finite 
amount of time. 

Thresholding ^ i ,u 

A corrected, standard deviation of the noise <r k for each subchannel k is 
estimated using a number of consecutively received symbols as described above in 
conjuntion with equations (10) and (11). A maximum number of bits to be allocated for 
each subchannel is then determined by comparing the estimated noise standard " r 
deviation cr k measured for each sxibchaiinel ^ to a series of r previously 

calculated and stored for each subchanneL Thus, rather than having to perform the 
complicated calculations in equations (8) and (9) including logarithms, exponentials, 
multiplications, etc., the present invention simply compares the estimated noise value 
measured for each subchannel to a number of thresholds. 

The thresholds themselves are calculated only once, (e.g., off-line at a 
factory), and are- subsequently used to establish a correspondence between a subchannel 
noise standard deviation value a k and a number . of allocated bits bfc for that subchannel. 
In other words, when the subchannel noise standard deviation value exceeds a certain 
limit, the number of bits that allocated to that subchannel changes. The subchannel 
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thresholds r are bbtained^by' solving equation (7) for a and inserting equation (8) 
which yields: ' "'"*'" "' '"' : ' tiJ 



*-(<?)= 



Is 



Q 



;s 2(2* M)- 1, : - 



(12) 



5 



where q . signifies. the discrete values at which the. subchannel signal constellation size 
changes,. In effect, the subchannel, constellation , size, changes with the measured 
subchannel noise^to maintain a, constant BEE,.- , The, thresholds define the values of 
measured subchannel. noise standard deviation ,cr .-, where. the constellation size changes 

- as* given by::.; ..v . .-; '.. ,.• " •.'...:.V:r.o.--.::; w. v .': - 

1 3 v 1 

10 ..where. ^ is the maximum number of bits to be .transmitted per subchannel, e.g., a 
maximum number of 12 bits. Once the thresholds are determined (preferably off-line) 
using equations (12) and (13) and the known values e„ 9 and Pm and the function Q(x) 
" : ■ for each subchannel, they are stored itf memory albrig with a corresponding number of 

" bits associated with each threshold: : ! ' " v ' r< ' - 

15 ;: If the measured .subchannel noise is so : strong,that not even .the smallest 

. ponstellation may be transmitted with the required bit error rate, that channel is not 
^ ,,,used, resulting in a zero bit allocation. Any subchannel with a bit allocation of zero 
reduces the number of active subchannels, by one. 

; The thresholding process is very simple compared to the complex 
20 calculations required, by conventional bit loading algorithms. The; thresholds are 
calculated following equations (12) and (13) and stored in memory. The received 
symbol mean error magnitude is measured for each : subchannel and used as an estimate 
of the subchannel noise standard deviation c k . The subchannel bit allocation is 
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associated to an interval between the two thresholds in which the value of the noise 
standard deviation estimate a k falls. The thresholding/comparison may be 
implemented in software using a rudimentary data processor, or preferably in relatively 
simply hardware to increase operational speed. 

5 ■ * ■ ' Fig. 7 is block diagram illustrating the thresholding operation. The 

measured subchannel noise standard deviation a k is input to a bank of parallel 

J - : comparators 70. Thd othef input of each coinparatof 70'is one of a number of 
■■ ' subbhanhel thresholds x r : ' calculated ta'accbfdahce with equation (12) for' different 
valties of 'g as indicated 4tte"4liatidn^i j)7 : Iii f m exarriple embodiment where b max is 12, 

mo : ' there are twelve comparators 70 and twelve 1 different stored 1 threshold -values' r, - r n 72. 
The combined outputs of each comparator 70 for a "bit vector" that is input to an. 
encoder 74. The least significant bit (LSB) of the bit vector is fixed to zero. In each 
comparator 70'; If the input noise deviation a k is less than the specific threshold value 
r , the comparator 70 generates a logical "1" at its output. Alternatively, if the input 

15 value is greater thairi the threshold, the comparator 70 generates a logical "0" at its 

{ ' /l output. ' ' :r: - l " 7L -•■ ••..?^:£i Jr,,M :<•"■ - : •<■-'"- 

As illustrated in Fig.. 8, t^he least significant bit in the comparator input bit 
vector supplied to the encoder 74 is in the zero, position, and the encoder ,74 outputs a 
zero if the bit vector contains only zeros. In effect, the specific comparator whose input 
20 threshold is first exceeded by the input noise standard deviation r M determines the 

"number of bills indicated by endbdbt ; 74 to : 'bfe allocated to that subchannel. In the case 
■ ! '- 1 shown in fig. 8, the smaller the noise* standard deviation value, the larger the number of 
bits that may be allocated to the subchartnel. " " 7 

Of course, the thresholding Operation 'illustrated in Figs. 7 and 8 is only an 
25 example. Another example* implementation is shown in Fig.' 9 where a single 
•xompafison device 84' is Used. The measured iioise standard deviation r k for a 
subchannel is 'stored' in an input register 80 and input id a comparator 84 . Digital 
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- threshold values r calculated for this subchannel are stored in threshold'memory 82 
and are sequentially input to the digital comparator 84 * The resulting sequential 
comparator outputs are passed through a demultiplexer in order to write different 
positions in the bit vector register 88. The encoder 90 analyzes the bit vector stored in 
register 88, and outputs the corresponding number of bits to be allocated for that 
subchannel based upon the thresholding comparison operation. 

If the total number of bits allocated to all subchannels in the thresholding 
step is less than the target number of bits, the desired bit rate cannot be achieved while 
imaintaining.the prescribed;bit error: rate -because- of increased noise or other factors. 
Consequently^ the data .transmission rate ; oflhe $ER performance- must be reduced. A 
!" decrease-in transmission rate :means^all6cating a fewer dumber of bits to each ] 
, subchannel. Alternatively, another set of thresholds r corresponding to a higher BER 
may. be employed. Thus, several sets .of threshold values for different BERs may be 
calculated using equations (12) and (13) and stored in memory. However, the former 
decreasejn transmission rate/allocated bits.is preferred.. ... 

Block Reduction 

; } If the total number of bits allocated across all subchannels, exceeds the 
target number of bits, excess,bits may be allocated, to the active subchannels such that 
the data transmission rate can exceed the r required data rate and still fulfill the 
prescribed bit error rate. If that goal is desired, the bit allocation procedures may end 
here. However, further bit allocation refinements may also be desirable. For example, 
the transmitted bit rate may be reduced to the prescribed bit rate to achieve certainty as 
well as 'an extra immunity to' noised 

If such further refinement is desired, a decision is. made if the total 
number of allocated bits exceeds the target number of bits by more than the current 
number of active subchannels, i.e., subchannels that have not been set to zero. If so, 
one bit is subtracted from the initially allocated number of bits to be transmitted on each 
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aetiye subchannel. Of course, any subchannels thathaye zero bit allocations; reduces 
the number of used subchannels ,by one: . ; . - ■ , . k •:, • 

This "block reduction" step is repeated until the number of allocated bits 
in excess of the overall target number of bits is less than the number of active 
subchannels k. The reduced number of allocated bits corresponds to a lower overall 
transmission rate with the margin evenly distributed among the subchannels. 



Fine Correction 



. \ r - A further fjne Qorr,eetida step.rnay be employed to compensate for the 

difference between the number .^f ibits remaiining to he allocated and/the target number 
of bus after the block reduction step, is 1 completed. JForeach active suhchannel^the 
:rckitii>n between the estimated noise standard deviation errand the threshold of its 
subchannel constellation is calculated 'as-: : r . : .c; -riT v--v" : ' - r ? ".tr.^ 



f 



where //ir[6 A .] signifies the operation of comparing the noise standard deviation a k of 
Subchannel k to the prbcalculaied thresholds for that subchannel. The distance dfc 
represents the difference between the" noise standard deviation value o k and the 
associated threshold r where the bitloading value changes. The subchannel having the 
largest distance is the first subchannel to have a bit removed from the groups of 
allocated bits. 

Because bits are subtracted from subchannels with the largest distance d k , 
the subchannel most sensitive to noise in its present constellation is advantageously 
reduced first. If the bit allocation of a subchannel that is not the most sensitive to noise 
is reduced, its immunity to noise is boosted. Unfortunately, the most sensitive 
subchannel remains unchanged thereby becoming the "weakest link". The distance 
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• - calculation ensures that- the weakest links are the first to have" bits removed then a 
. /relatively stroiig'subchannel.- - : ' '" ' ' •• - 

This task may take more or less time consuming depending on the total 
number of subchannels to be reduced'ahd on 'tlW choice of sorting algorithm. For 
5 v. example, assume there are 1024 subchannels khd that after the block reduction step 

, there are an excess of 1023 bits allocated, i.ei, 1023 bits more than the targeted number. 
l , : , • ..The task then is to reduce 1023 in arder of "largest distance first" by finding the 1023 
■.. largest elements in ^ the.' distance vector.? !^ on f the other hand, there is an excess of only 
10 bits, only the 10 largest elements in thb-distaric'e 'vector need be found which takes 
10 less time. The sorting algorithm which sorts the distance'vecitor elements from largest 
to smallest also affects this process. A fast sorting algorithm may be more complex to 
implement, require a sophisticated microprocessor, etc. 

A; flow chart diagram otitIMirig : exa^ 
, thresholding, block reduction, and fine cortectibh steps is how described' in conjunction 

15 wittUhe bit allocation -routine (block 1 00) allustl-ated ill Fig! 10: An initial series of 
operations is performed usually at the- time the device incorporating the present 
invention, such as a modem, is being manufactured: More specifibally, a series of 
threshold values r(q) is calculated for each subchannel in accordance with equations 12 
and 13 (block 102). In addition, each threshold value is correlated to a preset number of 

20 bits to be allocated for that subchannel. The threshold values and correlated numbers of 
bits are then stored. As mentioned above, different sets of threshold values . 
corresponding to bit allocations may be calculated for several different bit/symbol error 
rates. The calculations and procedures performed in block .1 02 are generally performed 
only once. 0 

25 In terms of determining the actual number of bits to be allocated for each 

subchannel, the noise deviation estimate a k is measured for each subchannel and 
corrected in accordance with equation (11) (block 104). In a preferred embodiment, the 



WO 99/16224 PCT/SE98/01598 

21 

actual correction factor is incorporated-in the one time threshold calculation process in 
block 102 since the correction factor is known. In,additioiji the number of active 
subchannels N sc is set to all of the subchannels (block 106). 

. The noise J standard devrntion estimate-for each subchannel . cr A . is 
compared to the. thresholds; r > predetermined for each subchannel from smallest to 
largest magnitude. When the)noise standard deviation first exceeds one of the 
thresholds, the stored, number of bits, /^corresponding to that first exceeded threshold is 
allocated to that subchannel, (block 108). For any subchannel where ho bits can be 
..allocated, the. subchannel is,deemsd r to, be inactive, and the number;'of active 
.subchannels..^ is decremented.: ?\vr<. f!v <;U:r«~. ■< j* ' 

A decisions made whether the sum of the bits allocated to all of the 
subchannels is greater than or equal to a target number of bits B target (block 110). If 
■ not > & e current channel^ conditions will not provide the targeted data bit rate at the 
.targeted bit error rate. A decision. must be. made then, to decide whethet to change the 
.bit error r rate requirements or the bititransmission rate, requirements (blocks - 
Depending upon .whether such ; reduced service requirements are acceptable, control 
may return to block 102^ or block 104. lV . - ., • r., t , , ,. . /. 

On the other hand, if the total number of bits allocated across all 
subchannels is greater than the target number' of bits," another decision is made whether 
the number of bits that can be allocated to the subchannels in excess of the target 
number of bits exceeds or is equal to the number of active subchannels (block 114). If 
so, one bit may be subtracted from the number of bits currently allocated to each active 
subchannel (block 1 16). Any subchannel whose group of allocated bits is reduced to 
zero after the subtraction operation causes a decrementing of the number of active 
subchannels N sc . This is an iterative process, with control returning to decision 
block 114 via the shown feedback path. ... « .. • • : 
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When the number of allocated bits in excess of the target number is less 
than the number of active subchannels N se , the distance dfc for each subchannel k is 
jie.term.ined in accordance with equation (J 4). (block 118). In order of decreasing , 
, distance.values for each subchannel k 9 one.bit is subtracted from each active subchannel 
5 ; yntil the sum, of the allocated bits, equals . the target number of bits (block 120). 

While the present method may be performed under relatively static' 
1 ' thannel conditions, it may aisb be advantageously erriployed in dynamic channel 
conditions. For example, the noise : 6fl : a^&fcnlaf subchannel may change over time. 
Accordingly, an inquiry is irirfde ih "ddcisiOii block 122 whether a subchannel condition 
|6 ' hks eKangcd. If ndt, c6ntinued ; period 

simply end. However, if the subchannel has changed, control returns to block 104 to 
obtain the updated noise deviation estimate for one or more subchannels affected by the 
condition change, and the method repeats?" ' ' - : 

The present invention may, of course, be applied to any 
1 5 1 telecommunications application that employes hiulti-tiarrief type modulation; e.g.,' DMT, 
"etc. "Thus, while! the invention has been 'described in connection With what is presently 
considered to be the most practical and preferred embodiment, jt is to be understood 
r that the invention is not to be limited to the disclosed -embodiment, but on the contrary, 
is intended to cover various modifications and equivalent r arrangements included within 

20 the spirit and scope of the appended claims. 

* .. ■ t , - t , . . . . . • . . ■ 
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WHAT IS CLAIMED IS: 

1 . In a multi-carrier communications system with a transmitter arid d 
receiver linked together by a commuriicatibris ctiaftnel affected by noise and other 
vaifying channel impairments where data bits to be transmitted over the channel kt a 
prescribed data transmission rate'arid bit errbr rate are'alloc&ted to plural' sutichahhels, a 
method comprising the steps, of . . ^ { , . 

calculating for each sufechannel a s.eriespf thresholds corresponding to a preset 
number of bits to be allocated to the subchannel; . . , . 

obtaining, an. error estimate pp,ejach ( sipbcliannfl;. .... , ?< , 

. comparing the error, -estimate for p^gh subchannel tp the thresholds r calculated for 
the. subchannel; and r ■ , ' = . , , 

allocating a number pf bits Jo each 
thresholds exceeded by the error estimate. ^ 

2. The methpd in claim 1 ? wherein the calculating step is performed only 
.., once and the obtaining, comparing, and allocating steps are performed for an initial 

condition of the channql and again when the initial condition of the channel changes. 

3 ! . "fh6 method in clkim 2/ furthef comprising: 

^erformirig the obtairiirig; comparing, and allocating steps in real time whbhever 
the channel condition changes. ; ' : : ' l 

4. The method in claim 1, wherein the error estimate is a noise standard 
deviation measured for each subchannel. 

5. The method in claim 4, wherein the noise standard deviation is corrected 
for nonuniform characteristics. 

6. The method in claim 1, wherein each subchannel is associated with a 
signal constellation where bits allocated to the subchannel are mapped to signal points 
in the signal constellation, and wherein the subchannel thresholds are calculated based 
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on the prescribed bit error rate, an average energy of the subchannel signal 
constellation, and the-srgnal -constellation size. > ' 



, ^ 7. r The method in claim 1 ? further comprising: 
. . determining whether a total number of bits allocated to all of the subchannels is 
greater than or equal to a total target number of bits associated with the prescribed data 
transmission rate. 



8. The method in claim 7, further comprising: 

if the total number of allocated bits is less than the total target number of bits, 
reducing the prescribed data transmission rate or reducing the prescribed bit error rate. 

9. The method in claim 7, further comprising: 

if the total number of allocated bits is greater than or equal to the total target 
number of bits, reducing the number of bits allocated to one or more of the subchannels 
until the total number of allocated bits is nearly equal to the total target number of bits. 

10. In a multi-carrier communications, system with a transmitter and a 
receiver linked together by a communications channel affected by noise and other 
varying channel impairments where data bits to be transmitted over the channel at a 
prescribed data transmission rate and bit error rate are allocated to plural subchannels, a 
method comprising the steps . of: : ;. , . ; \. : ■ ; 

obtaining an error estimate on each subchannel; - - : 

comparing the error estimate for each subchannel to a series of thresholds 
corresponding to a preset number of bits to be allocated to the subchannel; 

allocating a number of bits to each subchannel corresponding to one of the 
thresholds .exceeded by the^error estimate; and . 

when a total number of bits allocated to all of the subchannels is- in excess of a 
total target number of bits associated with the prescribed data transmission rate, 
reducing the number of bits to be allocated to each subchannel. 
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. 11. The method in claim. 10, further comprising:,: — : :. v.; 

determining an active number of subchannels having at least one allocated bit; 

and 

if the excess number of bits exceeds the active number of subchannels, reducing 
5 the number of bits to be allocated to each active subchannel by one bit until the excess 
number of bits is less than the number of active subchannels. 

12. The method in claim 1 1, further comprising: 

when a remainder of the excess number of bits is less than the number of active 
subchannels, further reducing the reminder by reducing the number of bits allocated to a 
ib selected one of the active subchannels. 

13. The method in claim 12, wherein the one active subchannel is selected 
based on a magnitude of a di fference between the error estimate and the exceeded.. 
threshold for each active subchannel, the active subchannel having the largest difference 
magnitude being selected as the one active subchannel. 

15 14. The method in claim 13, the further reducing step further comprising: 

subtracting one bit from each active subchannel in order of decreasing magnitude 
difference. 

15. The method in claim 10, repeating the steps in claim 10 when a condition 
of the channel changes. , / ^, 

20 16. The method in claim 15, wherein the steps in claim 10 are repeated in real 

time. 

17. The method in claim 10, wherein the error estimate is a noise standard 
deviation measured for each subchannel. - 

. . \, . •■ . » : l ' i A : ./ . V" v .. ' ' 8 * • 

18. The method in claim 17, wherein the noise standard deviation is corrected 

. • : * • • j ■ .' .ii « <:: "-ii ' . tic. .. ,4 .'! *•• • 
25 for nonuniform characteristics. 
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19. " The method inclaim 10, wherein each subchannel is associated with a 
signal c6nsteUation r wherebits allocated 'to the subchannel are encoded to signal points 
Withe signal constellation, and wherein the subchannel thresholds are predetermined 
based on the prescribed bit error rate, an average energy of the subchannel signal 

t constellation, and the signal constellation ?ize. 

20. A multi-carrier communications j sysftemxotnprisingl : 

a transmitter and a receiver linked together by a communications channel 
affected by noise and other varying channel impairments^ where data bits to be 
transmitted over the channel at a prescribed data transmission rate and bit error rate are 
allocated to plural subchannels; and 

electronic circuitry performing the following tasks: 

: comparing the error* estimate for each subchannel to a series of thresholds 
^corresponding to^presdtnumb v er of bits'to be allocated to the subchannel, and 
allocating a number of bits to each subchannel corresponding to one of the 
.:• thresholds exceeded by -the error estimate 1 .: t." .-i':; V:.j : f ;\ 

2 1 . The multi-carrier communication system according to claim 20, wherein 
the thresholds are calculated only once. and the comparing and allocating tasks are 
performed when the channel condition 'cfaihges:'-' • ; * r " : - 

22. The multi-carrier communication system according to claim 20, wherein 
the error estimate is a noise standard deviation measured for each subchannel. 

23 . The multi-carrier communication system according id claim 22, wherein 
the noise standard deviation is corrected for nonuniform characteristics^ ; 

24. The multi-carrier communication system according to claim 20, wherein 
each subchannel is associated with a signal constellation where bits allocated to the 
subchannel are entoded .to signal points in the signal -constellation, and wherein the 
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subchannel thresholds are calculated; based on th$ prescribed, bit error rate, an average 
energy of the subchannel signal constellation, and the signal constellation-size. 

25. The multi-carrier commiiriication system according to claim 20, wherein 
the electronic circuitry determines whether a total number of bits allocated to all of the 
subchannels is greater than or equal to a total target number of bits associated with the 
prescribed data transmission rate; t ./ r . • / ... :.. - 

26. The muiti-carrier communication system according to claim 25, wherein 
'if the total number of allocated bits is less tihan the total target number of bits, the 
electronic circuitry reduces the prescribed data transmission rate or reduces the 
prescribed bit error rate. 

. . .27. ; The multi-cairier communis 
ri if the : total 5 number of allocated t bitsas greater, than unequal tP the total, target number of 
, bits, the electronic: circuitry reduces thejnumber of bits allocated to. one or more of the 

subchannels until the total number of allocated, bits is nearly equaUo the total target 

number of bits. 

28, The rnukiTcarrier communication system according: to claim. 27, where the 
electronic circuitry performs the^jo^ving.ftuthte tasks:. 

determining an active number of subchannels having at least one allocated bit; 

and 

if the total number of allocated bits exceeds the total target number of bits by the 
. active number of subchannels or more, reducing-the number of bits to be allocated to 
each active subchannel rby one bit until the excess number of bits, is less than the number 
of active subchannels. 

, ; 29. < The multi-carrier communication system according to claim 28 r wherein 
wh^n a remainder of the excess number of hits is less than the number of active 



WO 99/16224 PCT/SE98/01598 

28 

subchannels, the electronic circuitry further. reducing, the remainder by reducing the 
number of bits allocated to a selected one of the active subchannels. 

30. The^ multi-carrier communication system according to claim 29, wherein 
the one active subchannel is selected based on a magnitude of a difference between the 
error estimate and the exceeded threshold for each active subchannel the active 
subchannel having the largest difference magnitude being selected as the one active 

subchannel/- ^ " \ :;; i y } 

~ ... ' ■ *• "• i 

I 31. The multi-carrier communication system according to claim 14, wherein 
the' electronic circuitry subtracts one bit from each adtiye: subchannel in order of 
decreasing magnitudc/differdnce. / t * n 

. 32. The multi-carrier communication system according to claim 20, wherein 
the electronic circuitry irfipleqients the tasks usingV processor. 

« • 33.. The..multircarrier communication system according to claim 20, wherein 

the electronic "circuitry implements the tasks using discrete hardware elements. 
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